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ABSTRACT. In this work, we investigated the effect of vaccina- 
tion on the transmission dynamics of rotavirus disease. Rotavirus 
can be spread both through direct or indirect transmission contact 
with infected children or unhygienic environment. Vaccination is 
administered in three doses. We determined both the disease free 
equilibrium point, Po, and the endemic equilibrium point, Py. 
The disease-free equilibrium point, Po, is both locally and globally 
asymptotically stable if the effective reproduction number Ry < 1 
and unstable if Ry > 1. The endemic equilibrium point, Přņy, for 
k = 0,1, 2,3 with vaccination exists if and only if Ry > 1. In case of 
no vaccination, the endemic equilibrium point, Poy, has a unique 
stable endemic equilibrium point whenever the basic reproduction 
number without vaccination Roy > 1 and a backward bifurcation 
exists. Numerical results show that vaccination reduces the degree 
of susceptibility and infectiousness when children are exposed to 
the disease. Hence help to reduce rotavirus disease among children. 


1. INTRODUCTION 


Rotavirus is the most common cause of severe gastroenteritis (di- 
arrhea) in infants and young children (Tate et al., 2012) leading to 
over 600,000 to 760,000 deaths annually worldwide (Clark et al., 2004; 
WHO, 2012; Gavi Alliance, 2013). Every infant is expected to be in- 
fected with rotavirus within the first years (0-5) of life (CDC, 2009). In 
the United States rotavirus infections affect approximately 2.7 million 
children under 5 years of age and result in the hospitalization of 55, 
000 children every year (Parashar et al., 1995). 
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Each year rotavirus causes an estimated 111 million episodes of diar- 
rhea requiring 2 million hospitalizations and 400,000 deaths in children 
under 5 years (Kim et al., 2010); 85% of these deaths occurring in low- 
income countries of Africa (Parashar et al., 2009). In Uganda, every 
year about 19,700 children under five die due to diarrhea caused by 
rotavirus (WHO, 2012). 


The primary mode of rotavirus transmission is fecal-oral (Hochwald 
et al., 1999). A confirmed case of rotavirus diarrhea can be defined 
in a patient with diarrhea containing rotavirus antigen detected by an 
enzyme immunoassay in a fecal specimen (Tate et al., 2012). Rotavirus 
can survive for months at room temperature (Romas et al., 1998) and 
can be passed from one person to another through a set of contam- 
inated hands with the virus or by touching a contaminated surface 
or object (Butz et al.,1993). The virus enters the body through the 
mouth. Children can spread rotavirus before and after they develop 
the symptoms (Parashar et al., 1995). 


Once infection occurs, the incubation period for rotavirus disease is 
about 1-3 days (Mastretta et al., 2002; CDC, 2009). Most primary 
rotavirus infections are associated with acute diarrhea that can lead 
to dehydration and occasionally to death. Common symptoms involve 
vomiting and diarrhea for 3-8 days, frequent fever and abdominal pain. 
Immunity after infection is incomplete but recurrent infections tend to 
be less severe (Bishop et al., 1973; Parashar et al., 1998). Gastroin- 
testinal symptoms generally resolve in 3-7 days (CDC, 2009). 


Many studies about rotavirus immunity have found that maternal an- 
tibodies protect younger infants (WHO, 2009) since breast milk con- 
tains lactadherine (Morrow et al., 2004) that reduces gastrointestinal 
infections. Lactadherine is thought to prevent symptomatic rotavirus 
infection (Newburg et al., 1998) when children are breast fed. Breast 
milk anti-bodies to rotavirus has however been cited as one of the pos- 
sible factors for the low rotavirus disease in the low income countries 
compared to the high income countries (Nelson et al., 1998) but does 
not kill the virus completely. 


Vaccination has proven to be efficacious in preventing rotavirus-related 
disease, gastrointestinal disease and health-care use in infants and tod- 
dlers (Vesikari et al., 2006; Vesikari et al., 2007; Ruiz-Palacios et al., 
2006). Two rotavirus vaccines are currently available: Rotarix and Ro- 
tateq (WHO, 2009a; CDC, 2013). Both vaccines are torelated, safe, 
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immunogenic and highly effective (Vesikari et al., 2004). They are ad- 
ministered orally from 2, 4, and 6 months of age (Pan, 2003). Other 
means like, improving hygiene, sanitation and access to clean water are 
not sufficiently effective in preventing rotavirus, (Butz et al., 1993). 


RotaTeq contains five antigens : G1, G2, G3, G4 and P1 (Kim et 
al., 2010). It is also a pentavalent bovine-human vaccine (Pan, 2003) 
that accounts for over 80% of rotavirus related diseases (Vesikari et 
al., 2004). The RotaTeq vaccine has been proved to have efficacy of 
approximately 75% against rotavirus gastroenteritis disease, and 100% 
against severe rotavirus related diseases (Molholland, 2004). 


Rotatrix is a single strain of human rotavirus (G1, P[8]) ( Pan, 2003) 
and it was licensed to GlaxoSmithkline (GSK) in Mexico in 2004. It 
is alive attenuated vaccine administered in 2 doses (WHO, 2009). The 
development of different rotavirus vaccines has been accepted globally 
as a number one priority to reduce children mortality as one of the 
Millennium development goals 2015 (Vesikari et al., 2004). 


The WHO has released a global recommendation that all countries 
include infant rotavirus vaccination in their national immunization pro- 
grams and the GAVI Alliance has promised to provide financial support 
for rotavirus vaccination programs to developing countries (Kim et al., 
2011; WHO, 2009). Severe dehydrating rotavirus infection occurs pri- 
marily among unvaccinated children aged 3-35 months (CDC, 2009). 
Further studies show that rotavirus vaccines are considered to be the 
most effective preventation method (Tate et al., 2012). 


Various mathematical models for the transmission dynamics have aided 
our understanding of the important factors driving epidemic patterns of 
rotavirus diarrhea in the community (Pitzer, 2009; Shim et al., 2001). 
These models can be used to estimate the expected direct and indirect 
effects of rotavirus vaccines (Pitzer, 2009; Van et al, 2010) and have 
generated predictions that agree with early observations of the impact 
of vaccination in the developed countries (Pitzer, 2009). 


1.1. Model Formulation. In this model, we have both vaccinated 
and unvaccinated children. The proportion of children that are unvac- 
cinated are denoted by (1— p) and these join the susceptible group S4 (t). 
Recruitment into this compartment is through birth by adults at a rate 
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A. Susceptible children from $;(t) when exposed, can acquire disease 
either from the different infected classes J,(t), 7 = 1,2,3 where Z(t) 
denotes infected children that are unvaccinated, J(t) denotes infected 
children vaccinated for the first dose, /3(t) denotes infected children 
vaccinated for the second dose or from environment class E(t). The 
force of infection of the unvaccinated children is defined by ~(S}, E, J;), 
where; 


3 


(1.1) (St, Bj) = $1 S0(6355) + 


j=l 


4,5,E 
K+E 





Basing on (1.1), €; measure the degree of susceptibility between Sı and 
I(t), vı measure the degree of susceptibility between Sı and E(t), 6; 
measures the degree of infectiousness in terms of probability contact 
with the different infectious classes, [;(t) and K denotes the level of 
pathogen concentration in the environment, E(t). Children in the in- 
fectious class [,(t) recover naturally at a rate a, while others die due 
to the disease at a rate dı. 


Furthermore we assume that vaccination is administered in three doses, 
that is, V1, denoting children vaccinated for the first dose, V2, denoting 
children vaccinated for the second dose and V3, denoting children vacci- 
nated for the third dose. The proportion of children vaccinated for the 
first dose is given by p. After a period of 4 weeks (CDC, 2009; WHO 
2009), children in V; either take the second dose or not. A proportion 
n takes the second dose at a rate 7, while (1 — 17) does not. Those that 
become susceptible join class S(t) at a rate kı. These children get 
exposed to infection at a rate ~(S2, E, Ij), where; 
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(1.2) Her B, I) = e282 O4 


j=l 


Va So E 
K+E 





From (1.2), €&2 measure the degree of susceptibility between S and J;(t), 
və measure the degree of susceptibility between So and E(t). Children 
in [j(t) die due to disease at a rate d2, while others recover naturally 
arate Qo. 


After the second dose, V2(t), children are vaccinated for the third dose 
or not. A fraction ¢ is vaccinated for the third dose at a rate Ta while 
(1 — @) is not. When the fraction (1 — ¢) is exposed they become sus- 
ceptible and join the S3(t) class at a rate kg. Children $3(t) become 
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infected at a rate w(S3, E, Ij), where; 





Infected children in J3(t) recover naturally at a rate a3 while others die 
due to the disease at a rate ds. 


We further note that the degree of susceptibility and infectiousness 
vary in different classes. The degree of susceptibility is measured by €; 
& vi, with i = 1,2,3 while the degree of infectiousness is measured by 
0j, j = 1,2,3. From the model description we take 6; > 02 > 63 and 
6, + 62+ 63 = 1. 


Children at V3(t) lose vaccine immunity at a rate 8 and they join 
the susceptible class S(t). All recovered children can lose immunity 
and join the susceptible group S,(t) at a rate e4. The pathogen pop- 
ulation is generated at a rate y while Infected children contribute to 
its enhancement through excretion at rates 01, 02, and o3 from I(t), 
I(t) and I3(t) respectively. The model flow diagram is shown in Figure 
(1.1). From the above assumptions, descriptions and the model flow 
chart together lead to the set of non linear- differential equations that 
describe the dynamics of the disease. 


dS, (t) 
dt 


dV, (t) 
dt 


dSo(t) 
dt 


dV2(t) 
dt 





=(1-p)A—4(S), E, L;) + eR — uSı + 8V3, 





= pA — (1 — n)k1 Vi — Nnm Vi — uV, 





= (1 a n)KiVi = (So, E, 15) — pS», 





= NT Vi — (1 — $)KkaV2 — T2 V2 — UVə, 
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FIGURE 1.1. A schematic of rotavirus disease with 
both unvaccinated and vaccinated children describing the 
mechanism of our model 
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where (5), E, Ii), Y(S2, E, I;) and W(S3, E, I) are defined in equations 
(1.1), (1.2) and (1.3) accordingly, with initial conditions: 

Sı (0) = Sio, vı (0) = Vo, S2(0) = So, V2(0) = Vo, V3(0) = 
V30, 53(0) = S30, L1(0) = tio, 
1,(0) = I, I3(0) = Izo, R(O) = Ro and N(0) = No. 


1.2. Invariant Region. Basing on the model (1.4), we assume that 
all the variables and parameters of the model are positive for all t > 0. 
The children population N(t) can be determined by 


dN (t) 
dt 
In absence of rotavirus, there is no death from rotavirus, that is, d = 0 
then 


= A — uN(t) -d(l + h + b) 


(1.5) —— <A puN(t) 
Solving for N (t) in (1.5), we get 


N(t) < ——[(A — uNo)e™ — Al, 
(1.6) CA oe 
> N(t)< 2 —( 

H H 





E 


A 
For No < —, from (1.6) we have 
u 


0 < N(t) < No 
so whatever the case N (t) is bounded above. 
That is, 
A 
(1.7) N(t) < N* = maz {No, =} 
u 


This implies that, if there is no disease, that is, J; = 0 for all 7 = 1, 2, 3, 
A 
N* = —, is the steady state population which is globally asymptotically 


stable. Similarly, considering the last differential equation in the system 
(1.4), that is, 


dE (t 
2. = VE — mE + 01, + cele + 0313, 
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and let the total pathogen population be E(t) , then we have 








(1.8) 

dE(t) 

BT a) E+ onl + ozh + osl < (y — m)E + (01 +0 + 03) N(¢), 
dE(t 

= EO < (y-m)E+ (or + or +os)NG) 

But from (1.7), we have N(t) < A which implies that 

(1.9) 

dE(t A 

a < (y= pn) E+ (01 +02 +03)Nce(t) < (y — m) E+ (01 +024 7s), 


Integrating both sides of (1.9) gives 


(0, +02 +03)A 
ulm — ) 





(1.10) E(t) < [1 + Be] 


where B is a constant. The pathogen population size E(t) becomes 


(0, +02 +03)A 
ulm — 7) 





0 < E(t) < 


provided ju, > ¥. 
Therefore, the feasible solution set of the pathogen population of the 
system (1.4) enters the region 





{Me = (E6) e RE < Gta 


ulm — 7) 


Therefore, the feasible set for our model system (1.4) is given by 


fo = (Sı, Vi, So, Va, 53, V3, 4, I2, Is, R, E) € R| 





(1.11) (S1, Vi, S2, V2, S3, Vz, hi, I2, I3, R, E) > 0; 
vos 4; By < ee 
ae ~ p(y) 


is a positively invariant set under the flow induced by the model system 
(1.4). Hence the system is biologically meaningful and mathematically 
well-posed in the domain 2. 


1.3. Model Analysis. We analyze and investigate both existence and 
stability of the equilibrium points, and the different reproduction num- 
bers. 
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1.3.1. Positivity of Solutions. 
Lemma 1. Let the initial data be 
{($1(0), Vi(0), $2(0), V2(0), S3(0), V3(0), Z1 (0), Z2(0), 13(0), R(0), E(0)) > 0} € R}. 
Then, the solution set 
{Si(t), Va (t), S2(t), V2(t), S3(t), V(t), H(t), b(t), Ist), RO), BO} 
of the system (1.4) is non-negative for all t > 0. 








Proof: 
From the first equation of the model system (1.4) 
dSı(t 
d3, (t) 
> —usS 
d < Mo} 


Integrating by separation of variables gives, 
dS, (t) / 
— > | —udt 
J S K 


S(t) Ss (Oe Srt > 0, 
and using initial conditions S1 (0) = S10, gives 


S(t) > Sioe ™ > 0, since u>0 


this gives 


Considering the second equation of model system (1.4) 








dV, (t 
ik Jee pA —(1—n)miVi — nV — pV, > —(1 — n)ki +91 + Ww) 
dV, (t 
ik ) > —((1— n)kı +7 + u) V 


Integrating on both sides by separation of variables gives, 


[HE > [--am tom + wae 





this gives 


Vı (t) > Vi(O)e J(G=n)«1+n7i+p)dt > 0 





since ((1—n)kı +N +u) >0 
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Similarly from the other equations of the system (1.4) we derive the 
following results: 


S(t) > Sa(0)e 7T "# > 0, since p> 0 





V(t) > Vo(O)eSC-d)natontHdt > Q since ((1—4)Ko+ bro +p) > 0 
V) > V(O S Ett > Q since (8B +p) > 0 
S3(t) > S3(0)e-J""* >0, since p>0 


tutd))dt > 0 since (ai +p+di)) >0 


iN 
pie 
— 
+ 
ee 
V 
MN 
= 
— 
oS 
Nat 
S | 
Se 
g 








N 





D(t) > L(0je7 (2 #+2))4# > 0, since (az + u+ d2)) > 0 














w 





L(t) > I(0)e7 J (3 t#+4))4 > 0, since (az + u+ ds)) > 0 


R(t) > RO)e ftt > 0 since (eg +p) > 0 


E(t) > E(0)efO-#)4# > 0 since (y — jn) >0 
Hence the components of the solution to system (1.4) are non-negative. 


1.3.2. Disease Free Equilibrium (DFE), Po. 

Let Po = (59, Vo, S2, Ve, S9, I}, I}, 19, R°, E°) be the disease free equi- 
librium point. This is obtained by setting the right-hand side of the 
model system (1.1) to zero. That is, 


dSi(t) _ dVi(t) _ dSo(t) _ dVa(t) _ V(t) _ dSs(t) _ R(t) 


dt dt dt dt dt dt dt 
This gives 


= 0 





0 = (1 — p)A — P(S1, E, Jj) — S51 + 8V3 

0 = pA — ((1 -nki +91 + u)v, 
(1.12) 0 = (1 = n)K1V1 — (So, E, Ij) — 1S, 

0 = nV — ((1— ġ)k2 + T2 + u) V2, 

0 = (1 — $)K2V2 — Y(S3, E, Ij) — 183, 
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with 4 = I = 13 = 0 & E = 0 and therefore, according to (1.1), (1.2) 
and (1.3), Y(S;, E, I;) = 0 for i, j = 1,2,3. 


Solving (1.12), the disease-free equilibrium point, Po, 
H = (SY, Vy, 53, Vy, Ve, Se, 0,0,0,0, 0) 




















where 
(1.13) 
go UPA Bor nT pn 
r u w(B+p) (1 — ¢$)k2 + Or + u) (1 -nri tam + u)’ 
Vo = pa 
(=i +71 +p)’ 
S9 = (1 = n)ki pA 
: u (=k tam +u) 
y? = NT på 
2 (= 6)ka +T + u) (Lm) +m +u)’ 
v? = PTa nT pA 
> (B+) (1 = d)K2 +T + u) (1 — nk +N + BH)’ 
90 = (1— d)ke NTI pA 





u ((l—$)ko + bt. + u) (1 — n)ki +971 +y) 


1.3.3. The Effective Reproduction Number with Vaccination, Ry. The 
effective reproduction number, Re refers to the average number of new 
infections generated by a typical infectious individual in a community 
where intervention strategies are in place. This is computed using the 
next generation operator approach as described by Van den Driessche 
and Watmough (2002) as follows. 

Let 


(i) F(x) be the rate of appearance of new infections in compart- 
ment 2. 
(ii) V;*(x) be the rate of transfer of individuals into compartment i 
by all other means. 
(iii) V; (x) be the transfer of individuals out of the compartment i. 


Then the disease transmission model consists of the system of equations 
x; = fix) = Fi(x) — Vi(z) 


where 
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The next important step is to obtain the disease-free equilibrium point 
Po. We then compute matrices F and V which are m x m matrices, 
where m represents the infected classes, defined by 





and 





OV; bet 

V= | P| with 1<i,j <m, 
Ox j 

respectively, F is nonnegative and V is a nonsingular M-matrix (a ma- 

trix with inverse, belonging to the class of positive matrices). Thus 

V~! is nonnegative and consequently FV~' is nonnegative. 


We then compute matrix FV~!, defined as the next generation ma- 
trix (Diekmann et al., 1990). The effective reproductive number with 
vaccination Ry is then defined as 


Ry = p(FV™) 


where p(A) is the spectral radius of matrix A (or the maximum mod- 
ulus of the eigenvalues of A). 

Applying the above method. From (1.4), after rearranging our set of 
equations from the infected class, we derive F; and V; as 





























€191 (AL, balə 0313) nEs, 
(1 14) bo €259(A11, AT 0313) wise 
; €353(011) + bola + 0313) + 2E 
0 
and 
(ay H di) l 
(az H dz) Ip 
1.15 Vi = 
( ) (a3 H ds) I3 











(u = VE — 011, — oala — 0313 


Obtaining the partial derivatives of (1.14) with respect to 11, Io, Is 
and E and evaluating at the disease free point Po gives 


€1510, €15102 €15103 an 

E€28201 €28202 €25203 ra$o 

€3:930, €35302 €38303 vis 
0 0 0 0 


(1.16) F= 
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Similarly, the Jacobian matrix of (1.15) works out to be: 
(1.17) 


(ay + +d) 0 0 0 

V= 0 (a2 + u + dy) 0 0 
a 0 0 (a3 + u + ds) 0 
O71 —02 =03 (m — 7) 


We then find the inverse of (1.17) which is used to compute FV~ and 
we have 

















(1.18) 
€15141 15101 615102 | S102 €15193 115103 vı Sı 
nı nıK(uı — 7) n2 noK ("1 — 7) n3 n3K(m —=7) K(m>7) 
€25261 25201 €25262 | 125202 €25203 128203 v2 S2 
FV =| m  mK(m-7) m  mK(m-7) m  n3K(m-7) Km- 
e3 5301 135301 €35392 | —_-¥3S302 €35393 135303 1383 
nı nıK (m — 7) n2 n2K(mı — 7) n3 n3K(m -=-7) K(m-7) 
0 0 0 0 
with 





nı = (a + u + di), 

Ng = (a2 + u + də), 

ns = (a3 + u + ds), 
at disease free equilibrium Py = (S9, V2, S9, V2, V2, S9, 0,0,0,0, 0). 
The effective reproduction number with vaccination, Ry, is computed 


as maz,{\; | FV} — AI |} = 0 for A; > 0,7 = 1,2,3. Thus, Ry, is 
given as 














(1.19) 
so vici so V202 so V303 
Ry = 2 (e 0i 4 ) pa (cot. H ) + —3 (cats H ) > ML >y 
V m YE Kur 1) m K(m — 7) n3 K(m — 7) i 


1.3.4. Analysis of the Effective Reproduction Number with Vaccination, 
Ry. From (1.19), let 














S! vioi S. V202 
Ti= 1 («0 H ) > Ye= 2 («202 + ) ; 
1 (ar Futd) O Kma) (a2 + w+) Km —7) 





s9 : 
T3 = 3 (cots t ae ) 
(a3 +u +d) K(m — 7) 


(a) Considering Yı we have 





0 
Yı = 1 (a0: + Le ) 
(ai +u +d) K(m-7) 


(1—p)A vici 
Ge + +d) (aa | K(m >) i 
( B PT2 NTI på (a0: : vici DE 
ular +u +d) (B+ p) (1 — 6)K2 + G72 + u) (0 - n)ki + 971 +y) K(m — 7) 
= Rov + Rav 
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where 
(1 — p)A ( V101 ) 
R €101 4 
ov ular + +d) | K(u1 -= 7) 
ae B PT2 nT1 pA (a0: pe Pidi ) 
(ai + utd) (B+ u) (1 — b)K2 + G72 +u) (1 -= n)ki +971 + 1) K(m — 9) 


We note that Roy is the basic reproduction number without 
vaccination and Rsy is the effective reproduction for children 
whose vaccine immunity wanes after the third dose and they 
become susceptible again. 

(b) Considering Tə 


S9 V202 
6 meg (0% K(m —7) 





E (1-—n)ki pA («ct P V209 ) 
ula + u +d) ((1— n)ki +71 + u) K(m — 7) 
= Riv 


representing the effective reproduction number of children who 
receive the first dose of vaccination but they are not taken back 
for the second dose. 


(c) Also considering Ts 





so ; 
T3 = 3 («cots | V393 ) 
(a3 + + d) K(u1- 7) 








_ (l-¢)k2 nT pA (cots 4303 ) 
Haztpta) (1 — 6)K2 + $72 + u) (1 —)K1 +T + y) ' K(u1 —7) 
= Rov 


We note that Ys denotes the effective reproduction number of 
children who receive the first dose and second dose but did not 
go back for the third. We denote Y; by Roy 


Thus the effective reproduction with vaccination Ry is given by 





Ry = Rov + Riv + Rov + Ray 
in which 


(i) Roy is the basic reproduction number without vaccination. 
(ii) Riy is the effective reproduction number for children who re- 
ceive the first dose of vaccination only. 
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(iii) Ray is the effective reproduction number for children who get 
first and second doses only. 

(iv) Rsv is the effective reproduction number for children who re- 
ceive all doses and when the vaccine waned, they joined the 
susceptible class again. 


1.4. Local Stability of the Disease-free Equilibrium Point with 
Vaccination, Pj. This can be determined by examining the linearised 
form of the system (1.4) at the disease free state Pp. The Jacobian ma- 
trix is computed by differentiating each equation in the system with 
respect to the state variables. we get 




















JP = 
(1.20) 
a 0 0 0 B 0 —e1 8101 —€15 02 —€15103 €4 =F uot 
0 =x 0 0 0 0 0 0 0 0 0 
0 xe -u o o0 0 —e28201 —e2S202 —e28S203 0 — 282 
0 NTI 0 =xX2 0 0 0 0 0 0 0 
0 0 0 T2 —x8 0 0 0 0 0 0 
0 0 o0 x7 0 pe 3 S30, —e3S302 —e3S303 0 2353 
0 0 0 0 0 0 X3 e1 8102 €1 8103 0 PEA 
0 0 0 0 0 0 €2$201 X4 €2:5203 0 v52 
0 0 0 0 0 0 €3:9301 €3:9302 X5 0 vas 
0 0 0 0 0 0 ai a2 a3 —x9 0 
0 0 0 0 0 0 oi o2 03 0 (y-m) 
with 


Xı = ((1—n)ki +07 + u), x2 = ((1 — ġ)k2 + O72 + u), 
x3 = & 8101 — (&ı + u +d), 


X4 = E2820 — (ag + u + d), X5 = €39363 — (a3 + u + d), X6 = 
(1 B n)Ka, XC (1 =T Q)k2, 


Xs = (8 + u), Xo = (u + €4) 


The local stability of P) is determined basing on the signs of the eigen- 
values of the Jacobian matrix (1.20). The disease-free equilibrium 
point, Po, is said to be locally asymptotically stable if the real parts 
of the eigenvalues are all negative, otherwise it is said to be unstable. 
For an eigenvalue à of matrix (1.20) we have | Jp, — Al |= 0. 


Clearly from (1.20), we can easily see that Ay = —p,A2 = —((1 — 
n)ki +771 + u), A3 = >H 

Ag = —((1-¢) kot b72+ 1), A5 = —(8+H), Ae = —y, and A7 = —(u+e4). 
The remaining four A's are evaluated from: 
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X3 — À  €1S182 €18103 a 
E29201 X4 =À €289203 v252 7 
(1.21) det NG, aes yo vib =0 
a Oy. Ga: tahini 


Computing for \’s of matrix (1.21), we get a characteristic equation of 
a fourth degree as: 








> 
A 
Qa 
= 
> 
w 
a 
N 
>` 
N 


(1.22) ee 


where 





ay =—(x3t [xa + (xs + (y — mı)))), 


V3 S' 
a2 = x3 (uas H1) — 93 s5) H 





v35.: v25: 
(xas ty — Ha) a >) €25203695900 — 022). 

















135: V2 5: 
a3 = — ( (xa + (Y = H1)) + (x5 (7 — Ha) — 03 s2) €28203€39302 — 02 232) ya- 
V3: 125 125 
(xas H1) — 93 = “= (e3830203 — 02X5) — €28203((7 — p1)e3 5302 — o2 2%) 








v35.: 125: 
a4 = (aosi — u1) — 93 =) + (€25263)(€3530203 — 02X5) — €25203((y — 11)€3S302 — o2 = )) X3 


Since we have a characteristic equation of degree four, we can ana- 
lyze this equation by employing the Routh-Hurwitz criteria for stabil- 
ity. Thus from (1.22) for Routh-Hurwitz criteria to hold, the following 
conditions have to be satisfied. 

(i) a, >0 

(ii) az > 0 

(iii) a4 > 0 

(iv) ayaga3 > a3 + aĉas 
The Routh-Hurwitz criteria is the necessary and sufficient condition 
for local stability of the system with all the eigenvalues have negative 
real part. Therefore if we analyze the above itemized conditions we 
will have, 


(i) a, > 0 if and only if (x3 + [x4 + (vs + (y — 1))]) is negative. 
(ii) a3 > 0 if and only if 

V3S3 
K 


12S: eos 
)4 ae (€3930203 — 02X5) — €252603((y — 141)€39302 — 2222) 








S 
)) €2:5203€35302 — 02 22 2) X= 


( (xa F (Y= #1) + (x5 (7 = 41) — 93 K 


v3 S3 
K 





(aosi H1) — 93 


it is negative. 
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(iii) a4 > 0 if and only if 


1353 


K 








128: 
(aoso — H1) -= 93 ) + (£29203) (€3830203 — oaxs)xa) > (e5280 — p1)E3 8302 — 02 = )) x3 


(iv) Since the above three critical conditions hold, after algebraic 
computations it can be shown that a,a2a3 > a + a?a4 holds. 


Since the first seven eigenvalues of matrix (1.20) are negative and also 
matrix (1.21) obeys the Routh-Hurwitz conditions shown above, thus 
the following theorem holds. 


Theorem 1.1. The disease-free equilibrium point Po of model system 
(1.4) is locally asymptotically stable if Ry < 1 and unstable if Ry > 1. 


1.5. The Global Stability of Disease Free Equilibrium Point 
with Vaccination (DFE), P. To determine the global stability of 
disease free equilibrium point with vaccination, P), we use the com- 
parison approach by (Diekmann et al., 1990), we consider the vari- 
ables representing the infected components and their rates of change 


of these variables. It follows that at disease free equilibrium point, Po, 
h = Í = Ts E 0, thus 














dl, (t 0 
| a vE (2 -= etn) 
h(t 

Ge) a ] E G=n)ki pA s9 
afto =(F-V) | iF a H (G=n)eitntitmyK KE 

3 3 0 

: (=¢)r nt pA s 

afto E(t) Yale ( uK = (=p terse) (G=n)K1+etaT1 +4) a 





ae 0 
dt 

where the matrices F and V are defined in equations (1.16) and (1.17) 

respectively. 


At disease free equilibrium point we have: 

















S? < UPA | B om nt AN 
KIBS p pn) O Atea n) (Caner F) for all t > 0. 
. . S2 (1— jk A 
Similarly, a < *—/™ (CREE Lee 
S$ (1—¢)K2 NTI pA 
Ke = (eet) aoa Or allt 2.0. 


Thus it can be seen that 











dl, (t) 
afte) h(t) 
(1.23) d fi i <(FHYV) ai 
ible Be 
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Since all the eigenvalues of the matrix (F — V) have negative real 
parts, it implies that (1.23) is stable if Roy < 1 and as t —+ co, we 
will have J; — 0, b —> 0, I — 0, and E&E —> 0. Therefore by 
comparison theorem in Lakshmikantham et al. (1989), it follows that 
(L, I, I3, Æ) —> (0,0, 0,0) and the three remaining equations of model 
system (1.1) gives us solutions as in (1.13). 

Thus (Si, Vi, So, Va, V3, S3, lh, Ip, Is, R, E) —> h ast — œ for Roy < 
1, implying that P) is globally asymptotically stable. Hence the fol- 
lowing theorem holds. 


Theorem 1.2. If Ry < 1, the disease free equilibrium point, P), of 
the model system (1.4) is globally asymptotically stable and unstable if 
Ry >1. 


1.6. The Endemic Equilibrium Point (EEP) with Vaccination, 
Py. The endemic equilibrium point with vaccination, Pi, k = 0, 1, 2, 3, 
where Poy, represents the EEP without vaccination, Piy, with first dose 

of vaccination, P3;,, second dose of vaccination and Py, third dose of 
vaccination, is obtained by setting the model system (1.4) to zero. This 

is done by expressing all our state variables in terms of the force of in- 
fection, that is, 


3 
(1.24) Y (S, E I) — (es: (0; I; ) + ae) rJ 1,2,3. 


j 
After algebraic computations, we get the following expressions: 
(1.25) 
* 1 * x * 
Si = pun ee aa + 8V5], 











A 1— 
Ve p ME e Canana 
((1 — n)ki +971 + u) Ll 
NTI PT2 
VS = Ve Vs = Vea 
. ((1 — ġ)k2 + O72 + u) ; : (6+ p) avo 


, 


u M 


Pa y e fo 
i (ai +u +d) 2 (a2 + u + d2) 2 (a3 + u + d3) 
alf + alz +asl3 ye _ (Oi + oal + o3]3) 
(ute) (um — 7) 








Hee 
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To establish the existence and stability of the endemic equilibrium point 
with vaccination, Pw, for k = 0, 1, 2, 3, of (1.4), we will examine four 
cases: 


) without vaccination, Pòy 

i) with first dose of vaccination, Piy 

) with second dose of vaccination, Py 
) with third dose of vaccination, Psy 


1.6.1. Case one: Endemic Equilibrium Point without Vaccination, Py. 
The endemic equilibrium point without vaccination, Poy, is obtained 
when %3 = %3 = 0. Thus (1.24) reduces to 





(1.26) wr = a S* (0 It) + oTr 
where 
1 E4 Qı 
S* = >A- y1 -M), M = , = 
i wh pil ~ (u+e€a4) (ai +u +d) 1 
1 Se i 
(a, + +d)’ (m — 9) 


Substituting S7, I¥, E* into (1.26) we get the following equation in 
terms of y7. 


(1.27) OIF WI) = Oi(Ad? + BY] +C) =0 


where 





i ee 
Ga yo+ +d)) \@ralliart pray Y 
Ke4a1 oid 





= | fi €4Q1 
2 ja a ' Gan tet) K) D e area i) 


8s (1- H] 
ulari +u +d) 


KNe161 void 
Baitp+d) n(m —y)(a1+n+4) 














uae G (a1 aa) eS (Cr a ae aa Cr is a) 


C can be reduced further to Roy — 1 as follows: 


A {O01 
a (argi) 
u(ay +p+d) K (“1 — 7) 
= Roy —1 
From (1.27), the endemic equilibrium point without vaccination, Py, 


has three roots, that is, when ~j = 0 giving the first root which corre- 
sponds to the disease free equilibrium point (DFE) and Ay}? + Buy + 
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C = 0 which has the remaining two roots corresponding to the endemic 
equilibrium point (EEP). These two roots are presented graphically as 
shown in Figure 1.2 and Figure 1.3. In Figure 1.2, we present the nature 


Force of infection versus RO for first root of the EEP 








unstable DFE stable DFE 


Force Of Infection 
è 
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0 0.5 1 0 0.5 1 
eps1 eps1 


FIGURE 1.2. When we plot force of infection (FOI) 
versus basic reproduction number, Roy (top plot), we 
note that as FOI increases also Roy increases, but when 
Ro = 0.25 the FOI becomes constant thus we except 
no more spread of the disease. It will only keep in an 
endemic steady state. Again, when we vary Roy ver- 
sus direct contact transmission rate (€,), these are di- 
rectly proportional to each other. Finally, we note that 
when we plot FOI versus c€; still FOI increases and when 
€, = 0.25, the disease goes into an endemic steady state. 


of the first root to have a unique equilibrium point with no backward 
bifurcation while in Figure 1.3, we present the nature of the second 
root which has both the stable disease free equilibrium (DFE) region 
and unstable endemic equilibrium point (EEP) region. We further note 
that Figure 1.3 has a backward bifurcation that occurs when the ba- 
sic reproduction number without vaccination, Roy = 1. This is clearly 
seen when we plot force of infection (FOI) versus (Rov). Therefore after 
the analysis of the two endemic equilibrium roots without vaccination 
the following theorem holds. 


Theorem 1.3. The endemic equilibrium point without vaccination, 


Poy, has a unique stable endemic equilibrium point if and only if Roy > 
1 
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x 10 Force of infection (FOI) versus RO for second root of the EEP 
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FIGURE 1.3. When we plot force of infection (FOI) 
against Rọ in the top plot, we note a stable disease free 
region when FOI= 0 and when Rọ = 1, the FOI starts 
to increase in the stable endemic region where we note 
that the disease starts to spread again hence a backward 
bifurcation. When Rp versus direct transmission con- 
tact rate (€,), these are directly proportional and when 
we plot FOI versus epsl, we still experience a backward 
bifurcation when ¢€, = 0.4. 


1.6.2. Case two: Endemic Equilibrium Point with First Dose of Vac- 
cination, Piy. The endemic equilibrium point with first dose of vacci- 
nation, Py, will be obtained when wy} = %3 = 0. This is because 7} is 
the force of transmission without vaccination. And w3 is the force of 


transmission under second dose of vaccination. 
Thus (1.24) is reduced to 


V295 E* 


(1.28) Wy = €293 (0213) + KLE 


substituting the expressions for 5%, I}, and E* from (1.28), the endemic 
equilibrium satisfy the following equation 


(1.29) PIF) = VAF + BY, + C) =0 
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where 


A= €29909 
B = naoz — (Mé2025502 — €202 KNno( u1 — Y) — UN2V202) 
g= (1—n)K1 ph 
= 
u (=n) tnm +u) 
C = unaKno(uı — y) — (He20253Kno(u — Y) + Unv28302) 
C can be reduced further to 1 — Rıy as follows: 


Gaja (23 i V20295 ) 
n Kno(p1 —7) 











= (l— nk, pA (« Gch. V209 
m = 242 meir a 

(as + p+ d)u ((1— n)ki +N +u) K(u1 — 7) 
= ] — Ry 


Solutions of (1.29) are v3 = 0 and f(43) = 0. wi = 0 corresponds to 
disease free equilibrium point (DFE) whose stability has been estab- 
lished under Section 1.3 and f (43) = 0 corresponds to a situation when 
the disease persists (endemic). In case of backward bifurcation, multi- 
ple endemic equilibrium must exist. This implies that equation (1.29) 
indicates that there are three cases we have to consider of f(w>) = 0 
depending on the signs of B and C since A is always positive. That is, 

(1) If B < 0 and C = 0 or B? — 4AC = 0 , then equation (1.29) 
has a unique endemic equilibrium point (one positive root) and 
no backward bifurcation possibility. 

(2) If C > 0, B > 0 and B? —4AC > 0, then equation (1.29) has 
two endemic equilibria (two positive roots), and therefore it’s 
possible for backward bifurcation to occur. 

However it’s important to note that C is always positive if Riy < 1 
and negative if Ryy > 1. 


Theorem 1.4. The rotavirus model with first dose of vaccination has, 
(i) Precisely one unique endemic equilibrium if C <0 = > Riy > 


(ii) Precisely two endemic equilibrium if C > 0, B < 0 and B? — 
4AC > 0 
(iii) None otherwise 


By this result, Theorem 1.4 gives a condition for existence of endemic 
equilibrium point with first dose of vaccination. 
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Theorem 1.5. The endemic equilibrium Py, with first dose of vacci- 
nation exists if and only if Riv > 1 


1.6.3. Case three: Endemic Equilibrium Point with Second Dose of 
Vaccination, P3),. The endemic equilibrium point with second dose of 
vaccination, Py will be obtained when ył = %3 = 0. This is because 
y} is the force of transmission without treatment and vaccination. And 
W93 is the force of transmission under treatment and first dose of vacci- 
nation. 

Again (1.24) reduces to: 

393 E* 

K + E* 

substituting the expressions for 5%, I}, and E* from (1.30), the endemic 
equilibrium satisfy the following equation 


(1.31) V3 (Ws) = ¥3(Avs’ + BY; + C) = 0 


where 





(1.30) $3 = €353 (8313) + 


A = €30303 
B = pngo3 — (163035303 — €303Kna( ui — Y) — N3V303) 
(1 — ġ)k2 NT pA 
u ((L— d)K2 + O72 +u) (A -= nei +N + y) 
C = uns Kns(uı — 7) — (wes0353 Ks (pa — 7) + wn3l25303) 
C can be reduced further to 1 — Roy as follows: 


C epics (25 | V30393 ) 
n3 Knz(m — 7) 


5? V303 ) 
=1-—— | e3034 
= (ebs K(m-7) 


Applying the same procedure used in analyzing equations (1.29) to 
(1.31), we note that C is always positive if Rəy < 1 and negative if 
Rəy > 1 leading to Theorem 1.6. 


S; = 











Theorem 1.6. The rotavirus model with second dose of vaccination 


has, 
(i) Precisely one unique endemic equilibrium if C <0 = Rəy > 
1 
(ii) Precisely two endemic equilibrium if C > 0, B < 0 and B? — 
4AC > 0 


(iii) None otherwise 
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Depending on this result, Theorem 1.6 gives a condition for existence 
of endemic equilibrium point with second dose of vaccination. 


Theorem 1.7. The endemic equilibrium Py, with second dose of vac- 
cination exists if and only if Rəy > 1 


1.6.4. Case four: Endemic Equilibrium Point with third dose of Vac- 
cination, Py. Having looked at the three cases above, in a similar 
way we can still verify the endemic equilibrium point for third dose 
of vaccination, R3y. Since the four reproduction numbers, that is, 


Rov, Riy, Rey and Ry come from Ry which we have already con- 
sidered. Thus 


Theorem 1.8. The endemic equilibrium Piy, with vaccination exists 
if and only if Ry > 1. 


2. NUMERICAL SIMULATIONS 


From our model (1.4) we verified some of the analytical results numer- 
ically . This was done using a set of parameters whose sources are 
mainly from Literature as well as assumptions. Some of the data was 
obtained from Uganda National private hospitals where rotavirus vac- 
cines are mostly given. We divided the model into four subsections, 
that is, 


(i) without vaccination (S1 RE), 
ii) with first dose of vaccination (V1 S212 RE), 
ii) with second dose of vaccination (V253/3RE), 
) we assume no or minimal infection after the third dose of vac- 
cination, 


such that we reduce the complexity of the model during simulation. 
We simulated the model using Matlab ODE solvers. However we only 
simulated (i) and (ii) since the results in (iii) showed minimal change 
from those in (ii). The numerical simulations of model (1.4) was carried 
out using the set of parameter values given in Table 2.1. 


2.1. Without vaccination (S;/4,RE). We assumed our initial con- 
ditions as $,(0) = 1000, /,(0) = 100, R(0) = 0, (0) = 100. We in- 
vestigated the effect of both direct and indirect contact rates on the 
transmission dynamics of rotavirus disease. Considering direct trans- 
mission we varied different rates of e} while for indirect transmission 
we varied vı. We further investigated the effect of shedding rate of 
infected children on the environment. 
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Table 2.1: Parameters values of the Model system (1.4) 
Parameter Description Value Source 
A birth rate 0.0018day~! 35] 
Pp proportion of unvaccinated 0.5 Estimated 
€ degree of susceptibility between Sı & I; 0.00052day—! Estimated 
m death rate 0.0018day—! 35] 
E2 degree of susceptibility between S2 & I; 0.00042day~ 1 5] 
€3 degree of susceptibility between S3 & I; 0.00032day~ 1 5,16] 
nı degree of susceptibility between S1 & E 0.002day~! 5,7] 
V2 degree of susceptibility between S2 & E 0.001day~! 5,14] 
V3 degree of susceptibility between S3 & E 0.0015day—! 8,9] 
a1 degree of infectiousness between S1 & I; 0.5 Assumed 
82 degree of infectiousness between S2 & I; 0.3 Assumed 
03 degree of infectiousness between S3 & I; 0.2 Assumed 
€4 natural immunity warning rate 0.0027day~ Í [6] 
B vaccine warning rate 0.00274day~ t! [13, 14] 
p rate of vaccination for dose one 0.060179day—! Estimated 
T1 rate of vaccination for dose two 0.0357day—! Estimated 
T2 rate of vaccination for dose three 0.0157 Estimated 
K1 rate of not receiving dose two 0.5 Assumed 
K2 rate of not receiving dose three 0.2 Assumed 
ay recovery rate from [4 0.2day~* 11,18] 
a2 recovery rate from T2 0.33day~t 11,18,32] 
a3 recovery rate from [3 0.667day—! 11,18] 
dı death rate without vaccination 0.0004466day~* 30] 
d2 death rate with first dose 0.00004466day~! Assumed 
d3 death rate with second dose 0.000004466day~! Assumed 
o shedding rate of I4 10 — 100cellsL—+ 1, 7] 
o2 shedding rate of I> 5 — 100cellsL—1 28,30] 
03 shedding rate of [3 0 — 100cells L71 1,8] 
Hı free-pathogen death rate 0.0667day~ 1 6] 
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FIGURE 2.1. Effect of direct transmission (€,) on the 
infected class 


2.1.1. Effect of direct transmission (€,) on the infected class. In Figure 
2.1 we note that an increase in the rate of direct contact rate (e1) 
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between susceptible children and infected ones, increases the number 
of infected children with rotavirus disease. Thus we should ensure that 
the level of contact of infected children is reduced as much as possible 
to lower disease transmission. 


2.1.2. Effect of indirect transmission (vı) on the infected class. In Fig- 
ure 2.2 it is noted that as we reduce the rate of indirect contact (1) 
between susceptible children and exposed environment, the number of 
infected children reduces while in case of high contact rate, (vı), the 
number of infected children increases. Thus we should fight to make 
our environments as clean as we can such that children do not con- 
tract the disease from dirty environments. This can be achieved by 
practicing proper hygiene and sanitation. 
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FIGURE 2.2. Effect of indirect transmission (v,) on the 
infected class 


2.1.3. Effect of the shedding rate (o1) on the environment class. In 
Figure 2.3 when the rate at which infected children shed off (c1) to the 
environment is low, there is less pathogens shed off to the environment 
and so less infection will be caused compared to when there is high 
shedding rate (c1) of infected child to the environment. 


2.2. With first dose vaccination (V,S2/.RE). Under this Section, 
we want to understand the effect of vaccination on the susceptible, 
infected and the environment classes. Here, we vary the rate at which 
children vaccinated for the first dose (%1) join these class. We assume 
that in case children are vaccinated, we expect many children to become 
less susceptible, few children in the infected class and low shedding rate 
to the environment class. 
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Effect of shedding rate on the environment 
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FIGURE 2.3. Effect of the shedding rate (o1) on the en- 
vironment class 


2.2.1. Effect of vaccination on the susceptible class. In Figure 2.4 when 
the rate at which, children vaccinated for the first dose leave the vacci- 
nated class to the susceptible class is little or small, very few children 
will move to the infected class and in case it is high, many are going to 
be exposed to the disease whom in the long run will become infected 
and join the infected class. Thus we should ensure to keep the rate 
(k1) as low as we can to reduce the degree of susceptibility to disease. 
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FIGURE 2.4. 
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Effect of kı on the susceptible class 


2.2.2. Effect of vaccination on the infected class. In Figure 2.5 if the 
rate at which vaccinated children move to the infected class is high, 
many children are going to become infected but if the rate is very low, 
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less children are going to be infected. This implies that, in case (K 1) is 
maintained as low as possible, the degree of infectiousness will be low 
and vice versa. 
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FIGURE 2.5. Effect of kı on the infected class 


2.2.3. Effect of vaccination on the environment class. In Figure 2.6 it 
is noted that if the rate at which children vaccinated for the first dose 
and joining the susceptible class Sy is very high, we expect many to 
become exposed to the disease who will later become infected. Due to 
infection, there will be high growth rate of pathogens in the environ- 
ment compared to when &, is very low. 
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FIGURE 2.6. Effect of kı on the environment class 
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2.3. Conclusions. This study formulated a mathematical model for 
the transmission dynamics of rotavirus disease. We included three 
different doses being administered as the case of Uganda. Both the 
disease free equilibrium (DFE) and endemic equilbrium points (EEP) 
were determined. Applying the Routh-Hurwitz criteria, the disease-free 
equilibrium point, Po, is locally asymptotically stable if, the effective 
reproduction number with vaccination (Ry), Ry < 1 and unstable if 
Ry > 1. The endemic equilibrium point, Pj, with vaccination ex- 
ists if and only if Ry > 1. Further more, in case of no vaccination, 
the endemic equilibrium point without vaccination, Py, has a unique 
stable endemic equilibrium point if and only if the basic reproduction 
number without vaccination, Roy > 1 and a backward bifurcation ex- 
ists as shown in Figure 1.3 when Roy = 1. Our numerical results 
show that both direct and indirect transmission contact rates play a 
high role in the transmission of rotavirus disease. Since vaccination 
helps to reduce the number of infected children thus it should be used 
to fight rotavirus among children. Again the environment should also 
be thoroughly cleaned to reduce the effect of children to environment 
transmission. 
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